x1=[2,-5,1];
v1=x1;
x2=[4,-1,2];

dot(x2,v1)
dot(v1,v1)

v2=x2-dot(x2,v1)/dot(v1,v1)*v1

norm(v1).^2
norm(v2).^2
